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The generation and manipulation of carrier spin polarization in semiconduc- 
tors solely by electric fields has garnered significant attention as both an in- 
teresting manifestation of spin-orbit physics as well as a valuable capability for 
potential spintronics device ^ | 2 | 3 | 4 l One realization of these spin-orbit phenom- 
ena, the spin Hall effect (SHE)^, has been studied as a means of all-electrical 
spin current generation and spin separation in both semiconductor and metallic 
systems. Previous measurements of the spin Hall effect™ 1 ™ 11 have focused 
on steady-state generation and time-averaged detection, without directly ad- 
dressing the accumulation dynamics on the timescale of the spin coherence 
time. Here, we demonstrate time-resolved measurement of the dynamics of 
spin accumulation generated by the extrinsic spin Hall effect in a doped GaAs 
semiconductor channel. Using electrically-pumped time-resolved Kerr rotation, 
we image the accumulation, precession, and decay dynamics near the channel 
boundary with spatial and temporal resolution and identify multiple evolution 
time constants. We model these processes using time-dependent diffusion anal- 
ysis utilizing both exact and numerical solution techniques and find that the 
underlying physical spin coherence time differs from the dynamical rates of spin 
accumulation and decay observed near the sample edges. 

Theories have predicte d 5 1 6 1 12 1 13 1 , and experiments confirmed^, that an electric current in 
a crystal with spin-orbit coupling gives rise to a transverse spin current via the spin Hall 
effect. Spin-dependent scattering of carriers by charged impurities (the extrinsic SHE)PEnU 
or the direct effect of spin-orbit coupling on the band structure (intrinsic SHE)P2H3] causes 
spin-dependent splitting in momentum space and a resulting pure spin current. Although 
not directly observable, the presence of this spin current can be inferred from the existence 
of non-equilibrium spin accumulation near sample boundaries. While extrinsic spin Hall 
currents generated by impurity scattering evolve on momentum scattering timescales (< 
1 ps), spin Hall accumulation is expected to develop on the much slower spin coherence 
timescale r (~ 1 ns). Since this timescale is of the same order as that desired for fast 
electrical manipulation of spin polarization in spintronics devices, understanding dynamics 
on this timescale is critical for both physical and practical insights to the extrinsic SHE 
processes. 

Steady-state observations of electrically-generated spin accumulate d 7 1 15 1 1 ^ are effective 
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for inferring r, but they cannot directly access the dynamical processes on the nanosec- 
ond timescale. In contrast, time-resolved spin dynamics with picosecond resolution are rou- 
tinely measured using ultra-fast optical pump-probe technique ^ 17 1 18 1 . Time-resolution of bulk 
current-induced spin polarization (CISP) was achieved using a photoconductive switcbP', 
but only precessional dynamics were observed due to short duration of the ultra-fast current 
pulse. Further, in contrast to SHE, CISP is a bulk phenomenon and consequently neither 
the steady-state^^ nor the time-resolvecP^ measurements investigated spatial dynamics near 
the sample edge. In the present work, we combine the spatial resolution afforded by scanning 
Kerr microscopjf^ with an optical probe pulse delayed relative to the electrical pump pulse 
to achieve both temporal and spatial resolution of spin polarization generated electrically 
by the SHE in an n-doped GaAs channel. Details of the experimental technique are shown 
in Figure 1 and are discussed in the Methods section below. 

Figure 2b-d shows the Kerr rotation^ (£) as a function of delay time t in a fixed magnetic 
field B applied along the y-axis. The spin polarization is generated by the SHE from a 
voltage pulse V(t) with length t p = 6 ns and amplitude Vb = 2 V (Fig. 2a). The laser is 
positioned at y = 126 fim so as to be close to the boundary spin generation with minimal 
clipping of the spot by the edge. During the pulse (0 < t < t p ), spin polarization builds 
up due to the SHE. After the pulse has passed (t > t p ) : the accumulated spin polarization 
undergoes decay and spin precession (Fig. 2b-d). We fit 0jc{t) for t > t p to an exponentially 
decaying cosine to extract an inhomogeneous depolarization time r* and Larmor precession 
frequency is L = gji B B/h where g is the Lande g-factor, \i B the Bohr magneton, and h is 
Planck's constant. Fits to vl(B) give \g\ = 0.346 ±0.002 (Fig. 2e), which is consistent with 
g expected for this doping leveP^. The depolarization time r* = 2.8 ± 0.1 ns measured at a 
specific spatial location should differ from the physical spin decoherence time r. Because spin 
polarization can diffuse due to accumulation gradients, there is a second pathway for spin 
depolarization beyond decoherence which depends strongly on the measurement location. 
We reconcile the dynamically measured r* with the intrinsic decoherence time r later in our 
discussion. 

Typical optical studies of electrically generated spin accumulation measure the time- 
averaged projection of spin precession as a function of applied transverse magnetic field B. 
In analogy with the Hanle effect of depolarization of luminescence, s z should depolarize for 
increasing B when 2i\v L ~ 1/r. 21 Therefore, coherence times in steady-state experiments are 
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typically extracted from the linewidths of s z (B). Near the sample edge, s z (B) is a Lorentzian 
lineshape analogous to the Hanle effectP, whereas it becomes more complicated away from 
the edges due to the interplay of spin precession and diffusioif^SD j n ^he current experiment, 
measurement of 9k(B) does not represent a time- averaged steady-state accumulation, but 
rather a snapshot at a fixed time t of the dynamic behavior of an electrically-generated spin 
ensemble in a magnetic field. 

Representative scans of 9k{B) at y = 126 /xm are shown for t — 1, 5, and 9 ns in Fig. 
2f-h with the magnetic field applied along the y axis. For small £, 9k(B) grows in a broad 
peak that narrows as t increases (Fig. 2f). Only for t <~ t p > r does 9k(B) approach the 
Lorentzian lineshape expected from a conventional Hanle analysis (Fig. 2g). For t > t p: 
6k(B) is primarily governed by spin precession, displaying characteristic periodic lobes of 
decreasing amplitude away from B = (Fig. 2h). 

In Fig. 3a, we use a longer pulse t p = 15 ns to investigate accumulation dynamics with 
the current flowing for various Vb- We fit 9k(B = 0) to an exponential saturation with a 
time constant r acc . For each Vb, r acc is around 40% of the r* measured from decay of the 
spin polarization (Fig. 3b). Both r* and r acc decrease weakly with Vb, which is expected 
due to electron heatingp^Zl. 

We characterize the magnetic field lineshapes by their inverse half- width B~h, which 
increases with t before quickly saturating (Fig. 3a). We can understand this evolution of 
B 1 / 2 in a simple physical picture^. Soon after the pulse turns on at t = 0, spins are all 
recently generated at the sample edge and have had little time for spin precession about 
B. For later times, spins have a larger spread in generation times (up to t p ) and have 
correspondingly more time for precession; hence, there is more depolarization for a given B 
and the Hanle curve narrows (Fig. 2g). For t >r, the average precession time is governed 
by r rather than t and the Hanle width becomes constant as in a steady-state measurement. 
The coherence times Ti/ 2 calculated from the saturation of B 1 / 2 near t ~ t p , agree with 
decay times r* and are consequently also longer than the accumulation times r acc (Fig. 3b). 
Diffusion analysis of the SHE accumulation is necessary to reconcile the observed differences 
between the timescales r acc , r*, and Ti/ 2 - 

Electrically-generated spin accumulation in GaAs can be modeled using drift-diffusion 
equation j 21 | 23 l We treat the channel as infinitely long since its length / is much larger than 
the spin diffusion length L s = 3.9 /xm found from steady-state measurements at B = 0. 
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This assumption reduces the problem to only one spatial dimension and precludes the need 
for general two-dimensional modeling including spin drift 11 . The SHE generates a spin 
current transverse to the in-plane electric field E = £x and proportional to the spin Hall 
conductivity (Tsh, j) — ^SH^ k Ek- The total current of the i spin component along y from 
both diffusion and SHE is f = —Dd y Si — asHE5 iz , where D = L 2 s /r is the spin diffusion 
constant. Spin conservation at sample edges y = ±w/2 is enforced with hard- wall boundary 
conditions normal to the edge (f y = —Dd y s l — asHE8 iz = 0). Diffusive boundary conditions 
accounting for spin-orbit effects at the sample edge (such as in the intrinsic SHE) would 
require modification for effects on the scale of the mean free patbpH. Assuming slow spin 
decoherence (relative to momentum scattering) at rate 1/r, the spin polarization will obey 
a continuity equation for &(y,t) including decay and precession terms: 

8^ dj i sHu t) 

^(y,t) = -fi(y,t) ~ + (W^)B x s( y ,t) (i) 

We first develop an intuitive picture of the time-dependent spin Hall processes using the 
exact solution to Eq. [j] under the simplest conditions. For £> = 0, the components of s 
are uncoupled and Eq. [j] can be solved using a Green's function to obtain an infinite series 
solution for s z (y,t). The diffusion equation for s z can be written as: 

Jti M ~ D W iV,t) + r M = FM (2) 

OE 

E(y,t) = V • (<7£7fE) = &sh-q^~ 

F(y, t) is a source function which contains all SHE terms. The homogeneous F = Green's 
function for Eq. [2] is: 

G(y, y', t, t') = Mk m y] sinMe^'^'eii - if) (3) 

m 

where m is an integer, 0(t) is the Heaviside step function, k m = (2m + 1)tt/w and A m = 
1/r + k^Ll/r. For time-independent integrating Eq. [3] yields a series representation 
of the one-dimensional steady-state solution to the spin Hall diffusion equation. 7 To obtain 
a time-dependent solution, we assume an ideal square electric field pulse of width t p and 
amplitude Eq, E = E [Q(y + w/2) — Q(y — w/2)][Q(t) — Q(t — t p )]. The corresponding source 
function is F{y ) t) = asHE[8{y + w/2) — S(y — w/2)][(Q(t) — B(t — t p )} and the solution is 
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found by integrating: 

/OO POO 
dy' dt'G(y,y',t,t') 
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T m (t) = <j (i_ e - A ^), 0<t <t p (4) 

e -A m (t-t p )^ _ e -A m tp^ t > t p 

The three regimes of the term- by-term time-dependence function T m (t) in Eq. [i]can each 
be observed in Fig. 2b. The lines in Fig. 3a, represent fits of 0x(t) to Eq. [I] keeping the first 
200 terms in Eq. [4] and convoluting the solution with the Gaussian profile of the laser spot. 
Since L s = 3.9 ± 0.2 /im was found independently from steady-state spatial measurements, 
the only fit parameters are r and an overall amplitude scaling. The best fit values for the 
parameter r are plotted in Fig. 3b and are significantly longer than the experimentally 
measured timescales r acc , r*, and Ti/ 2 . 

In Fig. 3c, we plot and calculations from Eq. [I] convoluted with the laser profile for 

y = 126, 124, 122, and 120 fim for Vb = 2 V using r = 4.2 ns obtained from the earlier fits. 
We fix the amplitude of the calculation from a fit to y = 126 /xm, and the remaining curves 
have no free parameters. For y away from the edge, 9 Kit) does not grow exponentially in 
t and we cannot define the time constant r acc as in Fig. 3a. Comparison of calculations 
from Eq. [4] with and without the spot size averaging reveals that the apparent asymmetry 
between the growth and decay times r acc and r* near the edge is primarily due to spatial 
averaging of these diffusion profiles over the Gaussian laser spot. The difference between r* 
and r is real, however; dynamically measured spin polarization near the sample edge evolves 
with a faster time constant than the underlying spin coherence time. 

We can understand the fast evolution of spin polarization from the interplay of diffusion 
and spin decoherence. Since polarization gradients cause spins to diffuse away from the 
sample boundary, spin depolarization must occur faster than decoherence of the electrically- 
generated spins. These dynamics are captured in the diffusion analysis of Eq. [4] by the fast 
decay rate A m of terms with large m in Eq. [4| Higher m terms are primarily responsible 
for the discrepancy between the best fit value for r and the faster timescales r acc and r* 
observed for spin accumulation and decay in Fig. 3b, but they only contribute significantly 
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to s z near y = w/2 where all terms are in phase. In this boundary region, timescales should 
differ most from the coherence time r. 

We numerically calculate ds z /dy at y = 126 /xm from spatial scans in Fig. 3d. The 
spatial derivative of s z (y) is proportional to the diffusive spin current and is non-zero at the 
sample edge in the presence of a compensating spin Hall current for < t < t p . After the 
spin Hall current disappears at t — t p , ds z /dy relaxes with time constant t - — 1 ± 0.1 ns to 
satisfy the diffusive j* = boundary condition (Fig. 3d, inset). Since V(t) evolves faster 
than r _1 , we expect the diffusive spin current at the sample edge to relax faster as well. 
Infinitesimally close to the boundary, ds z /dy should respond to changes in V(t) as fast as 
the momentum scattering time. Taking into account the finite laser spot size, we calculate 
Tj = 1.15 ns from our model, in agreement with the measured value. 

Introducing the magnetic field B = By couples the spin components s x and s z . For 
this regime, we perform numerical solutions to the system of coupled time-dependent linear 
differential equations represented by Eq. [T] using the best fit values for L s and r obtained 
from the earlier field independent analysis. For the numerical solutions, we use the exact 
pulse profile E(t) = V(t)/l measured from the oscilloscope (Fig. 2a) as a source. The curves 
in Fig. 2b-d and Fig. 2f-h are numerical calculations of and 6k{B). There are no free 

parameters except an overall scaling to match the amplitude of Ok- The full experimental 
data set 9(t, B) at y — 126 /im for t p = 6 ns is shown in Fig. 4a. Figure 4b shows the full 
calculation of spin accumulation from the numerical solution to Eq. [T]for y = 126 /im. 

The agreement between the experiments and calculations demonstrates that a single ho- 
mogeneous decoherence time r captures the various timescales observed in time-resolved 
measurements of the accumulation, decay, and diffusive dynamics of boundary spin polar- 
ization due to the extrinsic spin Hall effect. While diffusive timescales are set by the spin 
coherence time, evolution near sample boundaries can be limited by the faster response of 
the spin current. This spatial dependence of timescales could prove helpful for utilizing 
electrically-generated spin polarization in high-frequency semiconductor devices. 

Methods 

Channels of width w and length / are processed from a 2-/xm thick silicon-doped GaAs 
epilayer on 200 nm of undoped Alo.4Gao.6As grown on a semi-insulating (001) GaAs substrate 
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by molecular beam epitaxy (Fig. lb). The n-GaAs has doping density n = 1 x 10 17 cm -3 
and mobility \i = 3800 cm 2 /Vs at T = 30 K. The sample is mounted in a helium flow 
cryostat so that the channel (x-direction) is perpendicular to the externally applied in-plane 
magnetic field B (^/-direction) . A voltage V(t) applied across annealed Ni/Ge/Au/Ni/Au 
Ohmic contacts creates an in-plane electric field E(t) = V(t)/l along x. We desire an 
impedance of ~ 50 Vt in order to deliver the maximum broadband electrical power to the 
device; choosing w = 256 /xm and I = 130 fim yields a device with dc resistance R = 48 Q 
at T = 30 K. 

Time resolution of SHE accumulation is achieved by electrically-pumped Kerr rotation 
microscopy using a mode-locked Ti:sapphire laser tuned to 1.51 eV that emits a 76-MHz 
train of ~ 150-fs pulses. The pulse repetition rate is reduced to 38 MHz by pulse picking 
with an electro-optic modulator. Each laser pulse is divided into a trigger and a linearly 
polarized probe pulse. Spin polarization is generated at the sample edges by the SHE due 
to the current from a square electrical pulse of width t pi amplitude Vb, and 0.8 ns rise 
time applied to the sample from a pulse pattern generator triggered by the optical pump 
pulse. The linearly polarized probe beam is focused through a microscope objective to a 1 
/xm spot on the surface of the sample which can be scanned with submicron resolution. A 
balanced photodiode bridge measures the Kerr rotation of the linear polarization axis 6 k of 
the reflected beam which is proportional to the spin polarization along the z-axis s z . The 
leading edge of the electric pulse profile V(t) arrives at an electronically programmable delay 
time t before the arrival of the optical pulse. All reported measurements are taken at the 
center of the length of the channel {x — 0). 

An absorptive RF switch alternates the center conductor of the coaxial cable between 
the two complementary outputs of the pulse generator at frequency fy = 1.337 kHz. The 
RF switch only passes ac components of V{t) ) so the V baseline is restored by adding 
a square wave at frequency fy back onto the switched pulse train (Fig. lc) resulting in a 
modulation of the pulse amplitude +Vb and — Vb at frequency fy. 9 k is then measured with 
a lock-in amplifier analogous to the ac detection used in Refs. 1711101111 but with a definite 
phase relationship between electrical and optical pulses. 

The electrical pulse induces a time-dependent reflectivity modulation AR/Rof the optical 
beam during the pulse duration due to electron heatingpES that tracks the profile V(t) 
measured by an oscilloscope (Fig. 2a). We use this effect to calibrate t = and confirm that 
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the device acts as a proper 50 Q termination due to the minimal temporal pulse distortion 
at the sample. 
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FIG. 1: Experimental design for time-resolved measurement of electrically-generated 
spin polarization, a, Schematic diagram of time-resolved measurement of SHE (EOM, electro- 
optic modulator, BS, Beam splitter, PD, photodiode, FG, Function Generator), b, Optical micro- 
scope image of the sample with coordinate system defined (units in microns) showing the origin 
(black circle) and the location for most of the measurements (red circle). The bright yellow regions 
are the gold contacts and the grey region is the GaAs channel, c Illustration of the ac pulse scheme 
for lock-in detection. The pure ac components of the switch output (red) is added to a square wave 
at fy to create a triggered pulse train with amplitude modulation at fy w 1 kHz. 
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FIG. 2: Time-resolved measurement of the spin Hall accumulation, a, Voltage pulse 
profile V(t) measured by an oscilloscope (red line) and by reflectivity modulation AR/R (black 
diamonds). The arbitrary units of AR/R are scaled to match the scope voltage. The vertical 
dotted lines mark the time region < t < t p . b - d, Representative scans of Kerr rotation 9(t) for 
t p = 6 ns and y = 126 /im. The blue lines are calculations from Eq. [ljwith r =4.2 ns, L s = 3.9 
/im, and y = 126 /im. (e) vl extracted from cosinusoidal fits to Ox(t). f - h, Kerr rotation 9k(B) 
as a function of magnetic field at three representative times. The blue lines are calculations from 
Eq.0 
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FIG. 3: Dynamics of spin Hall accumulation, a, Kerr rotation 9k if) (upper panel), and 
inverse field width B~j 2 (lower panel) with t p = 15 ns at y = 126 fim for Vb = 2.0 V (black), 
1.5 V (red), 1.0 V (blue), and 0.5 V (green). The solid lines are fits from our model, b, The 
accumulation time r acc (green triangles), Hanle time Ti/ 2 (blue circles), decay time r* (black 
squares), and coherence time r (red diamonds) as function of voltage, c, 9k (t) for y = 126 fim 
(black), 124 fim (red), 122 fim (blue), and 120 /xm (green). The solid lines are calculations with 
amplitude fixed by a fit at y = 126 /im. d, Spatially-resolved Kerr rotation ^(y) near the edge 
for t = 5 ns (black), 7 ns (red), 9 ns (blue), and 11 ns (green). The inset shows the decay of 
jy oc ds z (y)/dy extracted at y = 126 fim after the voltage pulse is complete (t p = 6 ns). 
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FIG. 4: Modeling boundary accumulation from the spin Hall effect, a, Kerr rotation 
9x(B,t) at y = 126 fim for t p = 6 ns. (b) Calculations of s z {B,t) from Eq. [I] with r = 4.2 and 
L s — 3.9 /im at y = 126 /irn. 
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